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Abstract 

The short-distance expansion of the heavy- light currents qj^Q and 
q 7^75 Q is constructed to order 1/rriQ, and to next-to-leading order 
in renormalization-group improved perturbation theory. It is shown 
that the 10 x 10 anomalous dimension matrix, which describes the 
scale dependence of the dimension-four effective current operators in 
the heavy quark effective theory, is to a large extent determined by the 
equations of motion, heavy quark symmetry, and reparameterization 
invariance. The next-to-leading order expressions for the Wilson coef- 
ficients at order 1/mq depend on only five unknown two-loop anoma- 
lous dimensions, among them that of the chromo-magnetic operator. 
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1 Introduction 



Over the last few years, the heavy quark effective theory (HQET) has been 
established as a convenient tool to analyze the properties of hadrons con- 
taining a heavy quark |], 0, [|, f|, |5], It provides a systematic ex- 
pansion of hadronic matrix elements around the limit Aqcd/^q — > 0, in 
which one or more heavy quark masses tend to infinity and the effective 
low-energy Lagrangian of QCD exhibits a spin-flavor symmetry under the 
substitution of heavy quarks of different flavor or spin, but with the same 
velocity || [j| [L0[ [IT| . The existence of such a symmetry limit, and the estab- 
lishment of HQET to analyze the symmetry-breaking corrections to it, helped 
to remove much of the model dependence from the theoretical description of 
semileptonic weak decay processes of heavy mesons or baryons. 

The construction of a systematic 1/rriQ expansion of hadronic matrix el- 
ements consists of two steps: One first constructs the 1/rriQ expansion of the 
weak currents and of the effective low-energy Lagrangian using the operator 
product expansion and renormalization group techniques. This part of the 
calculation involves short- distance physics only and, in particular, is inde- 
pendent of the external states. In a second step, hadronic matrix elements of 
the HQET operators are parameterized by universal, mg-independent form 
factors. For the case of transitions between two heavy quarks, the short- 
distance expansion is known to order 1/rriQ, and the corresponding Wilson 
coefficients have been calculated to next-to-leading order in perturbation the- 
ory. The necessary hadronic matrix elements for meson and baryon decays 
have been investigated in detail up to order 1/mg. These results build a 
solid theoretical basis for a model-independent extraction of the Cabibbo- 
Kobayashi-Maskawa matrix element V c b from semileptonic B decays, with a 
theoretical uncertainty of only 5% (for a review, see Ref. |7[ and references 
therein) . 

The 1 / rriQ expansion is less developed for heavy-to-light transitions, which 
are, nevertheless, of great phenomenological importance: A thorough under- 
standing of exclusive decays such as B —>■ £ i? or B — > n £i? is necessary for a 
precise determination of V u b- For the most interesting processes of this type, 
the constraints imposed by heavy quark symmetry have been analyzed in the 



rriQ —> oo limit (see, e.g., Refs. |H| [TJ, [13], |14], |15[), but little is known about 



the symmetry-breaking corrections. In particular, already at order I/thq 
the short-distance expansion of the weak currents is only known in leading 
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logarithmic approximation |IEfl , and the hadronic matrix elements appear- 
ing at that order have only been analyzed for the simplest case, namely for 
heavy meson decay constants [ I7|| . In this paper we address the first problem 
and reconsider the short- distance expansion of heavy-light currents at order 
1 / rug beyond the leading logarithmic approximation. A detailed analysis of 
the hadronic matrix elements for B — > 7r £u decays at order I/tuq will be 
presented elsewhere 



After a brief review of the formalism of HQET relevant to our work, 
we calculate in Sec. § the Wilson coefficients appearing in the expansion 
of heavy-light currents at the one-loop level. This provides, in particular, 
the matching corrections between QCD and HQET at \i = tuq. In Sec. [|, 
we improve this calculation by using standard renormalization group tech- 
niques to sum the leading and next-to-leading logarithms [a s ln(m.Q / fi)] n and 
a s [a s ln(mg//i)] n to all orders in perturbation theory. This is a complicated 
problem in that one has to deal with a set of ten operators that mix under 
renormalization, with the technical complication that the 10 x 10 one-loop 
anomalous dimension matrix cannot be fully diagonalized. For the complete 
next-to-leading order solution of the renormalization group equation, one 
needs the anomalous dimension matrix at two-loop order. Although we do 
not calculate this matrix explicitly, we are able to infer much of its structure 
by analyzing the constraints imposed by various symmetries of the effective 
theory, as expressed in its Feynman rules, the equations of motion, and an in- 
variance under reparameterizations of the momentum operator. From these 
considerations, we can determine the 10 x 10 two-loop matrix up to only 
five unknown parameters. We then construct the exact next-to-leading order 
solution for the Wilson coefficients in terms of these parameters. In Sec. [5], 
we summarize our results and give a sample application. 



2 Short-distance Expansion 

In HQET, a heavy quark Q bound inside a hadron moving at velocity v 
is represented by a velocity-dependent field h v (x), which is related to the 
conventional spinor field Q(x) by || 

i + i> 

h v {x) = exp(im Q v-x) Q(x) . (1) 
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By means of the phase redefinition one removes the large part of the heavy 
quark momentum from the new field. When the total momentum is written 
as p = rriQ v + k, the field h v carries the residual momentum k, which results 
from soft interactions of the heavy quark with light degrees of freedom and 
is typically of order Aq CD . The operator |(1 + f) projects out the heavy 
quark (rather than antiquark) components of the spinor. The antiquark 
components are integrated out to obtain the effective Lagrangian |], 0, || |J || 

C cS = h v iv-D h v + \O kin + C mag (/i) mag l + 0(l/mV) , (2) 
2m Q L &J v 

where _D M = — ig s T a A^ is the gauge- covariant derivative. The operators 
appearing at order 1/rriQ are 

Okin = K {iD) 2 h v , mag = jh cr^G^K . (3) 

Here is the gluon field strength tensor defined by [iD 11 , iD u ] = ig s G^ u . 
In the hadron's rest frame, it is readily seen that Okm describes the kinetic 
energy resulting from the residual motion of the heavy quark, whereas O mag 
describes the chromo-magnetic coupling of the heavy quark spin to the gluon 
field. One can show that, to all orders in perturbation theory, the kinetic 
operator O k m is not renormalized Jl9|]. The renormalization factor C mag (/x) 
of the chromo-magnetic operator will be given later. 

In order to construct a systematic 1/rriQ expansion, one works with the 
eigenstates of the leading-order term in the effective Lagrangian and treats 
the 1 / rriQ corrections as a perturbation. The set of operators that appear at 
order 1/mg can be reduced by using the leading-order equation of motion 
iv-D h v = 0, since the physical matrix elements of operators which vanish by 
this equation are at least of order l/m 2 q. We have used this freedom to omit 
in (0) the operator h v {iv-D) 2 h v . 

The goal of this paper is the construction of the short- distance expansion 
for the heavy-light vector and axial vector currents = q^Q and = 
q 7^75 Q in terms of operators of the effective theory, beyond the leading 
order in 1/rriQ. We shall discuss the case of the vector current in detail. The 
general form of the expansion can be written as 

Ci(n) Ji + ^—H O j + J- £ AM T k + 0{l/m%) . (4) 
i 2m Q j Zm Q k 



3 



The symbol = is used to indicate that this is an equation that holds on 
the level of matrix elements. The operators {Jj} form a complete set of 
local dimension-three current operators with the same quantum numbers as 
the original vector current in the full theory. In HQET there are two such 
operators, namely 

J 1 = qYh VJ J 2 = qv^h v . (5) 

Similarly, {Oj} are a complete set of local dimension-four operators. It is 
convenient to use the background field method, which ensures that there is no 
mixing between gauge-invariant and gauge-dependent operators. Moreover, 
operators that vanish by the equation of motion are irrelevant at order I/tuq. 
It is thus sufficient to consider gauge-invariant operators that do not vanish 
by the equation of motion. A convenient basis of such operators is ||20| : 

1 = qj> t iph v , A = q(-iv- i D)'fh v , 

2 = qv^ipK, 5 = q(-iv-t))v»h v , (6) 
O z = qiD»h v , 6 = q(-tfr)h v . 

For simplicity we consider here the limit where the light quark is massless, 
m q = 0. Otherwise one would have to include two additional operators 
7 = m q Jy and 8 = m q J 2 . 

In (f|) we have also included nonlocal operators Tk which arise from an in- 
sertion of a 1/ itlq correction to the effective Lagrangian into matrix elements 
of the leading-order currents: 

T, = zJdyT{J 1 (0),O kin (y)}, 

T 2 = iJdyT{J 2 (0),O kin (y)}, 

T3 = z|d Z /T{J 1 (0),O mag (i/)}, 

T 4 = z|d 2 /T{J 2 (0),O mag ( 2 /)}. (7) 

The Wilson coefficients of these time-ordered products are simply the prod- 
ucts of the coefficients of their component operators, i.e. 

AM = CM, 

am = cm, 
am = cM c m&M , 

AM) = cMc ma M- (8) 



4 



Nevertheless, since these nonlocal operators can mix into the local operators 
Oj under renormalization (but not vice versa), it is convenient to include 



them as parts of the effective currents \W . 

Our goal is the calculation of the short- distance coefficients Bj(n) in (Q) 
beyond the leading logarithmic approximation. The coefficients Ci(fi) have 



been calculated at next-to- leading order in Ref. ||21|| . We shall perform the 
calculation using dimensional regularization with modified minimal subtrac- 
tion (MS) and anticommuting 75. This scheme has the advantage that, to all 
orders in 1/ttiq, the operator product expansion of the axial vector current 
can be simply obtained from (Hf) by replacing q — ► —575 in the HQET 



operators. The Wilson coefficients remain unchanged pOfl . The reason is 
that in any diagram the 75 from the current can be moved outside next to 
the light quark spinor. For m q = 0, this operation always leads to a minus 
sign. Hence it is sufficient to consider the case of the vector current in detail. 

Before we turn to explicit calculations, let us recall that there are non- 
trivial relations between some of the coefficients Bj(fi) and Cj(/i) imposed 
by a "hidden" symmetry of the effective theory, namely its invariance un- 
der reparameterizations of the heavy quark velocity and residual momentum 
which leave the total momentum unchanged. One can show that, as a con- 
sequence, operators with a covariant derivative acting on a heavy quark field 
must always appear in certain combinations with lower-dimension operators 



[18||. For the case at hand, there is a unique way in which the operators Oi, 
O2, O3 can be combined with Ji, J 2 in a reparameterization invariant form, 
namely 

V 2m Q / 2m Q 
f „ / ip \, O2 + 2O3 

This implies that, to all orders in perturbation theory, 

s x (/i) = cm , 

BM = 1 -bm = cm. (10) 

This is an important constraint, which has to be obeyed by any explicit 
calculation. 
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3 One-loop Matching 



The Wilson coefficients in (Q) are defined by requiring that matrix elements 
of the vector current in the full theory agree, to any order in 1/itlq, with 
matrix elements calculated in HQET. Order by order in perturbation theory, 
these coefficients can be computed from a comparison of matrix elements in 
the two theories. Since the effective theory is constructed to reproduce cor- 
rectly the low-energy behavior of the full theory, this "matching" procedure is 
independent of any long-distance physics such as infrared singularities, non- 
perturbative effects, and the properties of the external states. There is thus 
a freedom in the choice of the external states and the infrared regularization 
scheme, which can be exploited to simplify the calculations. We find it most 
convenient to perform the matching of QCD onto HQET using on-shell ex- 
ternal quark states and dimensional regularization for both the ultraviolet 
and infrared singularities encountered in the evaluation of loop diagrams.^] 
This scheme has the great advantage that all loop diagrams in the effective 
theory vanish, since there is no mass scale other than the renormalization 
point fi. This means that matrix elements in HQET are simply given by 
their tree-level expressions. We use momentum assignments such that the 
incoming heavy quark has momentum p = trig v + k (with 2v-k + k 2 /rriQ = 0), 
while the outgoing light quark carries momentum p' (with p' 2 = 0). Then, 
for instance, 

(Ji) = u g (p', sf)'fuh(y,s) , 

(0 1 )=u q {p'is') 1 ^u h {v 1 s), (11) 

etc., where u q (p' ', s') and Uhiv , s) are on-shell spinors for a massless quark q 
and a heavy quark in HQET, respectively. They satisfy]/' u q (p', s') = and 
$Uh(v,s) = Uh(v,s). To complete the matching calculation, one computes 
in the full theory the vector current matrix element between on-shell quark 
states at one-loop order. Taking into account that the relation between the 
heavy quark spinors in QCD and in the effective theory is H] 

u Q (p, s) = (l + u h (v, s) + 0{l/m%) , (12) 

1 The usefulness of this scheme was pointed out by Eichtcn and Hill in the original 
calculation of the renormalized effective Lagrangian at order 1 / raq Q . 
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we find at one-loop order 



2mQ/ 3n D — 3 



Jl + ^\_ (D _ 4) /, 2 + 5i±^ 



\ D-7 

X \ 2 \ _± ' 2m Q l y ~ ~'\~* ' 2m^ 

C-9J + 23 / Oj\ _ _ /a/v 2 / 0, 

D-5 \2m Q / v J \2m Q / \2m c 
, (13) 

where .D is the space-time dimension. The matrix elements on the right- 
hand side are evaluated in the effective theory. Note that the operators 
J\ and Ox, as well as J2, O2 and O3, appear precisely in the combinations 
required by reparameterization invariance, cf. @. Hence, our explicit one- 
loop calculation is in accordance with the general relations (fLOp. 

Next we expand the above expressions around D = 4 and subtract the 
poles in 

1 2 

- = j—p -7s + ln4vr, (14) 

corresponding to an MS renormalization of the operators in the effective the- 
ory, to obtain the following one-loop expressions for the Wilson coefficients: 

C 1 ( M )=5 1 (^) = l + ^(ln^-|), 
1 2a 

CM = bm = - bm = , 

B M = ^(sln^-l 

671 V fi 



3tt V /i 

BM = -p(m^-i). (15) 



37T \ /i 



The one- loop expressions for Ci(fi) have been first obtained in Ref. |21 |. Our 
expressions for Bj(fi) are new. 
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4 Renormalization Group Improvement 



4.1 Introduction 

For fi -C rriQ, the calculation presented above becomes unsatisfactory, since 
the scale in the running coupling constant cannot be determined at one-loop 
order. Yet a s {fi) and a s {mo) may differ substantially, and one would thus like 
to resolve the scale ambiguity problem by going beyond the leading order in 
perturbation theory. Furthermore, the perturbative expansion of the Wilson 
coefficients is known to contain large logarithms of the type [a s ln(m,Q / /j,)] n , 
which one should sum to all orders. Both goals can be achieved by using the 
renormalization group to improve the one-loop results derived in the previous 
section. 

Let us introduce a compact matrix notation where we collect the renor- 
malized operators Oj and T k , as well as the coefficient functions Bj(/j,) and 
v4fc(/i), into 10-component vectors 

6= (0 1 ,...,0 6 ,T 1 ,...,T 4 ), 

B = (B 1 ,...,B 6 ,A 1 ,...,A 4 ). (16) 
The Wilson coefficients obey a renormalization group equation 

(,iL_ f ) B - M = 0, (17) 

where 7* is the transposed 10 x 10 anomalous dimension matrix. It is defined 
in terms of the matrix Z that relates the "bare" operators to the renormalized 
ones: Ob are — Z O. In the MS scheme, Z obeys an expansion 

Z = l + J2^Z k (g s ), (18) 
k e 

with coefficients Z k (g s ) that depend on the renormalized coupling constant 
g s . Given this expansion, one can show that [^] 

d ~ 

*f=-9»-£- z i{j9s), ( 19 ) 

meaning that the anomalous dimension matrix can be computed from the 
1/e poles in Z. 
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The formal solution of the renormalization group equation reads 

B(ji) = Ufam Q )B{m Q ), (20) 
with the evolution matrix |23], |24| 

U(fx,m Q )=T g ex V J &9j^jr- (21) 
9s{m Q ) 

The /^-function f3(g a ) = /idg s /dfj, describes the scale-dependence of the renor- 
malized coupling constant in QCD. The symbol "T 5 " means an ordering in 
the coupling constant such that the couplings increase from right to left (for 
< mq). This is necessary since, in general, the anomalous dimension ma- 
trices at different values of g do not commute: [7(^1), 7(^2)] 7^ 0. Eq. (pl|) 
can be solved perturbatively by expanding the /3-function and the anomalous 
dimension matrix in the renormalized coupling constant: 

%>)= ^+^(t0 2 + -... (22) 

Here 

A) = H-|n/, /^102-yn/, (23) 

where n / denotes the number of light quark flavors with mass below wiq . At 
next-to-leading order in renormalization-group improved perturbation the- 
ory, the evolution matrix can be written in the form |24] 



Ufa mo) = {l - U fa m Q )[l + ^^s) + ..., (24) 



where 

^m^exjj-l.^} (25) 

describes the evolution in leading logarithmic approximation, and S con- 
tains the next-to-leading corrections. This matrix is defined by the algebraic 
equation 

* + 2s^=2k*-w*- (26) 



9 



Finally, we expand the initial values of the coefficient functions at // = ttlq 
as 

B(m Q )=B + ^±B 1 + .... (27) 



From our one-loop results in ( |I"5| ) we obtain 

B = (1,0,0,0,0,0,1,0,1,0) 



Bx = \ (-8, 4, 8, -8, 24, 8, -8, 4, 5, 4) , (28) 



where we have used that [0 



26 a B (m Q ) 
3 4vr 

Putting everything together, we obtain from (|2Q| ) 



a iag (m Q ) = l + ^^^. (29) 



£(//) = Ufa, m Q ) [b + (fii + 5^0 

-^S'^o(^m Q )S + .... (30) 

What remains to be calculated are the matrices Uq and S. To this end, 
we need to know the 10 x 10 anomalous dimension matrix 7 at two- loop 
order. Note that the leading-log evolution matrix, as well as the coefficient 
(Bi + S Bo) that multiplies a s (mg), are independent of the renormalization 
scheme ||23|| . The matrix S that multiplies a s (/i), on the other hand, is 
scheme-dependent. Only when the /i-dependent terms in B(fi) are combined 
with the //-dependent matrix elements of the renormalized HQET operators 
does one obtain a renormalization-scheme invariant result. 



4.2 Anomalous dimension matrix 

Let us now analyze the structure of the anomalous dimension matrix in more 
detail. We will see that, to a large extent, the texture of this matrix can be 
determined without an explicit calculation. From @ and ( |10D we know that 
the operators 0\, O2, O3 and renormalize multiplicatively. They can mix 
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into 04, O5, Oq, but not vice versa. Therefore, the anomalous dimension 
matrix must be of the form 




7 = 7 B I , (31) 

with diagonal matrices 



7h i = diag( 7 hl ,7 hl ,7 hl 



7D = diag( 7 hl , 7 hl , 7 hl + 7 mag , 7 hl + 7 mag ) . (32) 

Here 7 hl is the universal anomalous dimension of dimension-three heavy- 
light current operators in HQET, and 7 mag is the anomalous dimension of 
the chromo- magnetic operator in (§). The form of ju is a consequence of 
reparameterization invariance [cf. ([HJ)], whereas the structure of 70 follows 
from (§). The one-loop coefficients of the anomalous dimensions are JlO, || 



7o hl = -4 , 7o mas = 6 . (33) 

The two-loop coefficient 7] 11 = — ^ — p7r 2 + Y n S nas been calculated in 
Refs. ipij, p5fl , whereas 7™ ag is not yet known. 

The 3x3 submatrix 75 in ([J]]) may be constructed by noting that the 
equation of motion iv-D h v = can be used to rewrite 

4 = — iv ■ d Ji , 

05 = — iv ■ d J2 , 

0,-0, = -id»[qh v ]. (34) 

The total derivatives of the currents on the right-hand side renormalize mul- 
tiplicatively and in the same way as the dimension-three operators Jj. The 
additional power of the external momentum carried by the current does not 
affect the divergences of loop diagrams. It follows that 

(7*)^ = 7% + M7aV (35) 

Let us then focus on the 3x3 submatrix ^a, which describes the mixing 
of Oi, 2 , 3 into 4 , 5 , Oq. To compute X4, one has to calculate in the 



11 



effective theory the ultraviolet divergent 1/e poles of the p'-dependent terms 
in the matrix element of the generic operator qYiD a h v , where p' denotes the 
external momentum of the light quark. At the end of the calculation, one 
substitutes T = 7^7", v^j 01 , g^ a for Oi, 2 , 3 , respectively. The structure 
of these terms is strongly constrained by the form of the Feynman rules of the 
effective theory, and by the equations of motion. In HQET, a heavy quark 
couples to gluons proportional to its velocity |3], ^ , and hence there do not 
appear any Dirac matrices on the right-hand side of T. Furthermore, since 
we consider the limit of a massless light quark, there can only be an even 
number of 7-matrices on the left-hand side of T, since vertices and propaga- 
tors come in pairs. Moreover, the equation of motion for the heavy quark, 
iv-D h v = 0, requires that the matrix elements must vanish upon contraction 
with v a . Likewise, the equation of motion for the light quark, ilf) q = 0, im- 
plies that the p'-dependent terms must vanish when one substitutes T = 7 Q r'. 
Using these constraints, we find that, to all orders in perturbation theory, 
the ultraviolet divergent pole in the matrix element must be of the form 



' q T iD a h v ) 



pole 



3e 



U q (p') 3p' a - 2v -p'v a - V -p'^ai) Y U h {v) 



+ p'-independent terms, 



(36) 



where the function fi(g s ) has a perturbative expansion in powers of the 
coupling constant. From this, it is straightforward to derive that 



7a = 7 



2 

' 3 


1 

3 



■1 

"3 


2 
" 3 



2\ 



1/ 



(37) 



where 



dh 
dg s 



(38) 



We can use similar arguments to impose some restrictions on the 4x3 
submatrix 7^, which describes the mixing of the nonlocal operators into 
O4, O5, Oq. The Feynman rules of HQET imply that 



X 2 > 



f2(g s 



pole 



v ■ p u q (p') Y Uhiy) + p'-independent terms, 
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+ ^'-independent terms, 



(39) 



where T = 7 M for 7\ and T 3 , whereas r = v M for T 2 and T 4 . It follows that 



Id 



7 b 
7 C -47 c - 27 

Vo -3 7 







(40) 



J 



where 



t 8/2 
1 = ~ 9 'W. 



~9s ~7\ j 7 = ~9s - 



r= - g -w/ r= - g -w/ (41) 

We have calculated the anomalous dimension matrix 7 at the one-loop 
level by computing the ultraviolet divergences of the bare operators Ob are - 
Our results are in agreement with the general relations derived above. Ex- 
panding the anomalous dimensions 7°, . . . , 7 d as in (22), we find the one- loop 
coefficients 



7o=4 



7o 



7o = 7o 



(42) 



32 

~~3~' ' u ' u 3 
We note that the first row in ja, as well as the first and third rows in 7c, 
can be derived from a previous one-loop calculation by Falk and Grinstein 
[ jn| . We confirm their results. 



4.3 Leading-log solution 

Given the one-loop anomalous dimension matrix, one can calculate the evolu- 
tion matrix in leading logarithmic approximation from (fZ5p. The evaluation 
of this equation would be straightforward if there were a matrix V such that 
V~ 1a /qV were diagonal [24fl. However, such a matrix does not exist in the 
present case. The best one can achieve is to construct a matrix W which 
brings into Jordan form: 



W- 1 ? q W = ij. 



(43) 



13 



A matrix which accomplishes this is: 



W 



( 





























1 \ 






















1 





-1 




























1 










i 











4 


32 


16 


1 


64 


2 




3 


27 


3 


27 


3 


27 


3 




2 


16 


32 


16 


8 


128 





26 


32 


52 




3 


9 


3 


9 


9 


9 


3 


27 


3 




1 











4 

3 








-1 





-2 






















1 






















1 








1 

3 

























1 

















V 





1 








4 
3 














) 



(44) 



A Jordan matrix is convenient enough for an exponentiation in closed form. 
In our case, we only need that 



exp 



a b 
a 



exp (a) 



1 b 
1 



It is then straightforward to compute the evolution matrix from 



U (fi,m Q ) = W exp 



7 J ln <M«iQj 



(45) 



(46) 



,2/3 Q "~ a 8 (fi) 

Multiplying Uq with the tree-level matching condition as encoded in B in 



»D, we recover the leading- log results of Refs. |]16 

2/A) 



CM 

Bs(ji) 



x 



BM 
B 2 ((i) =BM = 



34 x 2/(3 ° 




4 x~ 1/(3 ° 




10 


16 


27 




27 




y- 






+ 


88 x- 1/Po 




20 




27 


27 




y 




~2x 2 ^ - 


4 


x -i/p + 


10 








3 


~3~ 







X 



X . 



with 



x 



«s(/-t) 

ocs{m Q ) 



(47) 
(48) 
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For mQ pa /i, these expressions can be expanded using 

(49) 

and one readily recovers the logarithmic terms in the one-loop coefficients 
given in ([15]). 

At this point, it is worthwhile to compare the the leading logarithmic 
approximation with the one-loop results. For the purpose of illustration, 
we use a s {mo) = 0.20 and a s (/i) = 0.36 (corresponding to mq pa m b and 
\i pa 1 GeV), as well as n/ = 4. From fl47D , we then obtain: = 1.15, 
= 5 ll = Qj 5 ll = o.63, £LL = _ .38, E^l = -0.21. In the one-loop 

expressions (^), we use the average value a s = 0.28, as well as toq/jU = 4.8. 
This gives: = 1.02, = 0.06, 5^ = 0.12, b[ 1) = 0.44, = -0.02, 
5^ = -0.07. Obviously, there are significant differences between these two 
approximation schemes, which can only be resolved by going beyond the 
leading order. 



4.4 Next-to-leading order solution 

Let us then discuss the construction of the complete next-to-leading order 
solution of the renormalization group equation. According to (|30|) , we need 
to construct the matrix S that satisfies the algebraic equation (p6|). To this 
end, it is convenient to define T = W SW, where W has been given in 
(H). The algebraic equation which determines T reads 



T + ±r [jj, f] = ^- W~^\ W-^jj. (50) 

We have solved this equation be modifying the Jordan matrix, 7j — > jj(rj), 
such that 7j(0) = Jj, but Jj(t]) can be diagonalized for 77 7^ 0. It is then 
straightforward to construct a matrix T(rj) which satisfies the above equation 
with 7j replaced by 7,7(77) ||24|| . The desired matrix T is obtained by taking 
the limit 77 — > 0. We find that most of the entries in T vanish. Among the 
nonvanishing components are: 

Tu = Shi 4" S mag j i = 2, 5, 

7;* = Sk; i = 3,4,6,7,8,9,10, (51) 
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where 

These are nothing but the next-to-leading corrections to the Wilson co- 
efficients of the dimension-three heavy-light currents and of the chromo- 
magnetic operator, which are given by |^T], 0, ||: 

C Ol) = x 2/ Po L + a a (m q ) - a a (/,) _ 4 a^l 

I 47T 3 7T J 

3 7T 

cw/0 - ^ 3/g ° (i + °- (mq ' " a,W *U + 1 • < 53 > 

L 47T 6 7T J 

As previously, x = a s (/i)/a 3 (mg). The remaining nonvanishing components 
of T are: 



T 15 



7i+7i 

4 7l mag - 4 7l " + 3 7l d 
6(/5 - 1) 



T 9 16 7l mag + 27 7l c + 9 7 f 

T 32 - - i T 35 - -3T 65 - ig2(/3o _ 3) , 

T34 = T67 = -64^-2l- (54) 



Given these results, one can compute the Wilson coefficients from (|30|) . 
For the coefficients of the nonlocal operators T^, and for the first three 

of the coefficients Bj(n), we confirm the general relations (Q) and fllQl) at next- 
to-leading order. The exact next-to-leading order results for the remaining 
coefficients are rather complicated. We present them in the form of three 
independent combinations: 



45 4 + B 5 + 2B 6 = ™ Cl l-L lux - T 34 ^l^iM) 
3 1 2/3 4tt J 

2 f 24 1 

+ -C 2 |l + 8C mag + — lnxj, 
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64 

3B 5 + 2B & = — — C\ 

y 



32 



-a 



mag 



mag 



+ 



47T 

29 3 , 

2 Ci — 2 C2 — - Ci C mag + - T 15 — 
10 '1 



1 + 62] 



32 



4tt 



+ fc 2 k 

4 



+ T 3 



15 



Arc 



as(m Q ) a s (m Q ) - a a (ji) 

+ 1 li — 



3tt 



47T 



(55) 



where we have omitted the /i-dependence of the coefficient functions for sim- 
plicity. The coefficients Cj(/x) have been given in 



5 Summary and an Application 

We have presented a detailed analysis of the short-distance expansion of the 
heavy- light currents q 'j^Q and q 7^75 Q to order 1 / vtlq in the heavy quark 
expansion, and to next-to-leading order in renormalizat ion-group improved 
perturbation theory. Our main result is that the 10 x 10 anomalous dimen- 
sion matrix 7, which describes the scale dependence of the dimension-four 
effective current operators in the heavy quark effective theory, can be deter- 
mined to a large extent from symmetry considerations. By evaluating the 
constraints that arise from reparameterization invariance, the equations of 
motion, and heavy quark symmetry, we have shown that, to all orders in 
perturbation theory, 7 can be expressed in terms of the universal anomalous 
dimension of the leading order (dimension-three) currents, the anomalous 
dimension of the chromo- magnetic operator, and four functions ^ l {g s ) of the 
coupling constant. We have calculated these functions to one- loop order, and 
have treated their two-loop coefficients, as well as the yet unknown two-loop 
anomalous dimension of the chromo-magnetic operator, as free parameters in 
the exact next-to-leading order solution of the renormalization group equa- 
tion. Our final expressions for the Wilson coefficients appearing at order 
l/m,Q in the expansion of the currents are given in (g), (|TUD, and ([55]). 

We believe that the fact that our results depend on only five unknown 
two-loop anomalous dimensions (as compared to the a priori 100 unknown 
entries in the two-loop anomalous dimension matrix) makes it feasible to 
obtain the complete next-to-leading order solution for the Wilson coefficients 
in the near future. This would be desirable since, as we have pointed out, 
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there are substantial numerical differences in the results obtained in leading 
logarithmic approximation as compared to those obtained from a one-loop 
calculation, which we have presented in Sec. |3|. These differences can only 
be resolved at next-to-leading order. 

The short-distance expansion of heavy-light currents considered in this 
paper plays an important role in applications of the heavy quark effective 
theory to decay processes such as B — > iv or B — > X£v, where X is a 
light meson (X = tt, p, etc.). As an example, we briefly discuss the case of 
meson decay constants, following the analysis of Ref. [pL7| ■ To order I/tuq in 
the heavy quark expansion, the decay constants of pseudoscalar and vector 
mesons can be written as 



4 



x 1 + 



1 

m Q 



A 



GM--b(p) 
6 



2d M 



(56) 



+ 



m Q 



A 



Cmag(/i) G 2 (m) - J^ B (V) 



where &m = 3 if M is a pseudoscalar meson, and du = —1 if M is a vector 
meson. F(/j,), Gj(/i), and A are mg-independent low-energy parameters of 
the effective theory. The coefficients b(n) and B(p) contain combinations of 
the Wilson coefficients Bj(n), which we have calculated in this paper. One 
obtains [O 



B 



KCi + C 2 ) 
3 

4(Ci + C 2 ) 



B\ — Bi + B4 + i?5 + B§ 



B\ — B 2 + -B4 + B 5 + S 6 



4d 
1 



B\ — B 3 — 3-B 4 — B e 



(57) 



B\ — B 3 — 3B 4 — B t 



Using the next-to-leading order expressions for the Wilson coefficients given 
in (1551), we find that 



16 k g a (fi) + a,(mq) _ g « s (m Q ) -a a (//) 



(3 a s (m Q ) 



71 



71 



(58) 



7 
9 



41 8 \ a s (m Q ) 
27 3 3 V tt 
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From the fact that the physical decay constants fu in ®) must be \i- 
independent, one can deduce the scale dependence of the low-energy pa- 
rameters. At leading order in l/rnq, we recover the well-known result that 
Ci([i) F(n) is \i- independent. At order 1/mg, the scale independent combi- 
nations are — |A6(yu) and C maig (fj,) G 2 (/i) — ^AB(/j). Note that the 
structure of the coefficients b(fi) and B(fi) is consistent with this require- 
ment, i.e., the /^-dependent terms in b(/j,) are independent of itlq, and those 
in B(fi) are proportional to C mag (/i). This provides a nontrivial test of our 
results.0 As a consequence, we can define renormalization-group invariant, 
fi- and mg-independent low-energy parameters by 



GT = GM - - \—]na a (ji) + 8T 3 , 



01 



Ml 



7T 



(59) 



G 



ren 
2 



a 8 (ji) 



-3/A) 



4A 



1 + ^T 32 ^ 

2 7T 



In terms of these parameters, one obtains, for instance, 

GM -±b(ji) = GT + I [j o \na s (m Q ) - (l - 8T 34 ) , (60) 

and a similar relation for the invariant combination involving G%(il). These 
relations show explicitly the nonanalytic mg-dependence of the I/tuq cor- 
rections to meson decay constants. 
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2 We have checked that this "consistency" of our short-distance analysis with the re- 
quirement of renormalizability of hadronic matrix elements is also obeyed for B — > tt £ v 
transitions, as well as for the decay constants of excited mesons. 
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